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We present a study of the decays — >■ DK^ with D mesons reconstructed in the K~^tt~tv'^ or 
K~tv^tt'^ final states, where D indicates a D" or a I?" meson. Using a sample of 474 milhon BB 
pairs collected with the BABAR detector at the PEP-II asymmetric-energy e'*"e~ collider at SLAC, 

we measure the ratios = p|g±I^[j^I^T^o]° ■ We obtain R^ = (5j;j^Q(stat)]'^^(syst)) x 10^'^ 
and R' = (l2iig(stat)+5(syst)) x W'^ , from which we extract the upper limits at 90% probability: 
R^ < 23 X 10~^ and R~ < 29 x 10^'^. Using these measurements, we obtain an upper limit for the 
ratio rs of the magnitudes of the b ^ u and b ^ c amplitudes rs < 0.13 at 90% probability. 



PACS numbers: 13.25.Hw, 14.40. Nd 



I. INTRODUCTION 



CP violation effects are described in the Standard 
Model (SM) of elementary particles with a single phase 
in the Cabibbo-Kobayashi-Maskawa (CKM) quark mix- 
ing matrix Vij [1]. One of the unitarity conditions 
for this matrix can be interpreted as a triangle in the 
plane of Wolfenstein parameters [2], where one of the 
angles is 7 = arg{-l/Jj,Krf/V^hV;d}. Various methods 
to determine 7 using — > DK^ decays have been 
proposed [3~5]. In this paper, we consider the decay 
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channel B+ DK+ with D K-tt+tt" [6] studied 
through the Atwood-Dunietz-Soni (ADS) method [4]. In 
this method the final state under consideration can be 
reached through 6 — c and h u processes as indicated 
in Fig. 1 that are followed by either Cabibbo-favored or 
Cabibbo suppressed D° decays. The interplay between 
different decay channels leads to a possibility to extract 
the angle 7 alongside with other parameters for these 
decays. 

Following the ADS method, we search for 
[/!r~7r+7r°]£ii4r+ events, where the favored 5+ D^K+ 
decay, followed by the doubly- Cabibbo-suppresscd D'^ 
A'^TT+TT^ decay, interferes with the suppressed — > 
D^K^ decay, followed by the Cabibbo-favored D° — 
K^TT^'K^ decay. These are called "opposite-sign" events 
because the two kaons in the final state have oppo- 
site charges. We also reconstruct a larger sample of 
"same-sign" events, which mainly arise from the favored 
D'^K^ decays followed by the Cabibbo-favored 
I?" _R'+7r^7r° decays. We define / = K^7T~Tr^ and 
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FIG. 1: Feynman diagrams for — s- D (top, h 
transition) and D^K^ (bottom, b —>-u transition) 



f = K tt+tt'^. We extract 



R- 



T{B+^[f]DK+) 
r(B + ^[/]BA' + )' 



r(i 



T{B-^[f]DK-)- 



(1) 
(2) 



from the selected and B~ samples, respectively. 

While our previous analysis [7] used another set of ob- 
servablcs: 



^ADS 
^ADS 



^ r(B + ^[f]DK+)+T{B-^[f]oK-) 
T(B+^[f]DK+)+T{B--^[f]DK-)' 

^ T{B- ^[f]DK-)-T{B + ^[f]DK+) 
~ T{B + ^[f]oK+)+T(B-^[f]oK-) 



(3) 
(4) 



we prefer to use observables defined in Eqs. 1 and 2 since 
their statistical uncertainties, which dominate in the final 
error of this measurement, are uncorrelated. 

The amplitude of the two-body B decay can be written 



A{B+ D°K+) = \A{B+ ^ D"K+)\rBe 

\AiB+^D"K + )\ 



(5) 



where tr = L ,„ — is the ratio of the masnitudcs 

\A{B+^D°K+)\ ^ 

of the b ^ u and b ^ c amplitudes, (5b is the CP conserv- 
ing strong phase, and 7 is the CP violating weak phase. 
For the three-body D decay we use similarly defined vari- 
ables: 



r {D" ^ f) 



J dm Alc^irn) 



^ J rfmADcs(m)AcF(m)e'^(") 
^Jdm Al^g (to) / dm A^p (m) 



(6) 
(7) 



where Acp(m) and ADcs("i) are the magnitude of the 
Cabibbo-favored and doubly-Cabibbo-Suppressed ampli- 
tudes, respectively, 6(m) is the relative strong phase, and 
TO indicates a position in the D Dalitz plot of squared in- 
variant masses [w^^,m^^o]- The parameter ku, called 
the coherence factor, can take values in the interval [0, 1]. 



Neglecting Z?-mixing effects, which in the SM give neg- 
ligible corrections to 7 and do not affect the measure- 
ment, the ratios i?+ and i?~ are related to the B- and 
_D-mesons' decay parameters through the following rela- 
tions: 

R+ =r%+r% + 2rBrDkDCOs{-f + S), (8) 
R- =r%+rl+2rBrDkDCoa{-f~S), (9) 

with 5 = 6b + Sb- The values oi ko and Sd measured 
by the CLEO-c collaboration [8], ko = 0.84 ± 0.07 and 
6d ~ (47^17)°, used in the signal yield estimation 
and tb extraction. The ratio has been measured in 
different experiments and we take the average value = 
(2.2 ± 0.1) X 10"'^ [9]. Its value is small compared to 
the present determination of r^, which is taken to be 
(0.106 ± 0.016) [10]. According to Eqs. 8 and 9, this 
implies that the measurements of ratios R^ is mainly 
sensitive to tb- For the same reason, the sensitivity to 7 
is reduced, and therefore the main aim of this analysis is 
to measure i?"*", and r^. The current high precision 
on is based on several earlier analyses by the BaBAR [7, 
11-13], BELLE [14-16], and CDF [17] collaborations. 

This paper is an update of our previous analysis [7] 
based on 226 x 10^ BB pairs and resulting in a measure- 
ment of i?ADS = (13l}o) ^ 10~^ wl^ich was translated 
into the 95% confidence level limit rs < 0.19. 

The results presented in this paper arc obtained with 
431 fb"^ of data collected at the T{AS) resonance with 
the BABm detector at the PEP-II e+e" collider at SLAC, 
corresponding to 474 x 10^ BB pairs. An additional "off- 
resonance" data sample of 45 fb~^, collected at a center- 
of-mass (CM) energy 40 MeV below the T{AS) resonance, 
is used to study backgrounds from "continuum" events, 
e+e~ qq {q = u, d, s, or c). 



II. EVENT RECONSTRUCTION AND 
SELECTION 

The BaBAR detector is described in detail else- 
where [18]. Charged-particle tracking is performed by a 
five-layer silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH). In addition to providing precise position 
information for tracking, the SVT and DCH measure the 
specific ionization, which is used for identification of low- 
momentum charged particles. At higher momenta pions 
and kaons are distinguished by Cherenkov radiation de- 
tected in a ring- imaging device (DIRC). The positions 
and energies of photons are measured with an electro- 
magnetic calorimeter (EMC) consisting of 6580 thallium- 
doped Csl crystals. These systems arc mounted inside a 
1.5 T solenoidal superconducting magnet. Muons are 
identified by the instrumented flux return, which is lo- 
cated outside the magnet. 

The event selection is based on studies of off-resonance 
data and Monte Carlo (MC) simulations of continuum 
and e+e- T(iS) BB events. The BaBAR detec- 
tor response is modeled with Geant4 [19]. We also use 
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EvtGen [20] to model the kinematics of B meson de- 
cays and JetSet [21] to model continuum background 
processes. All selection criteria are optimized by max- 
imizing the S/y/S + B ratio, where S and B are the 
expected numbers of the opposite-sign signal and back- 
ground events, respectively. In the optimization we as- 
sume an opposite-sign branching fraction of 4 x 10"^ [9]. 

The charged kaon and pion identification criteria are 
based on a likelihood technique. These criteria are typ- 
ically 85% efficient, depending on the momentum and 
polar angle, with misidentification rates at the 2% level. 
The 7r° candidates are reconstructed from pairs of pho- 
ton candidates with an invariant mass in the interval 
[119,146] MeV/c^ and with total energy greater than 
200 MeV. Each photon should have energy greater than 
70 MeV. 

The neutral D meson candidates are reconstructed 
from a charged kaon, a charged pion, and a neutral pion. 
The correlation between the tails in the distribution of 
the isTTTTr" invariant mass, and the 7r° candidate 

mass, m^o , is taken into account by requiring Imjj — m.^o | 
to be within 24 MeV/ of its nominal value [9] , which is 
1.5 times the experimental resolution. 

The B~^ candidates are reconstructed by combining D 
and candidates, and constraining them to originate 
from a common vertex. The probability distribution of 
the cosine of the B polar angle with respect to the beam 
axis in the CM frame, cos 0b, is expected to be propor- 
tional to (1 — cos^ Ob). We require | cos 6*5 1 < 0.8. 

We measure two almost independent kinematic vari- 
a bles: the beam-energy s ubstituted mass mEs = 
\/[s/2 +po ■PbY/E'^ — Pb^, and the energy difference 
AE = Eb ^ ■v/*/^, where E and p are the energy and 
momentum, the subscripts B and refer to the candi- 
date B meson and e^e~ system, respectively, -y/s is the 
center-of-mass energy, and Eb is measured in the CM 
frame. For correctly reconstructed B mesons the distri- 
bution of Toes peaks at the B mass, and the distribution 
of AE peaks at zero. The B candidates are required to 
have AE in the range [-23, 23] MeV (±1.3 standard de- 
viations) . We consider only events with toes in the range 
[5.20,5.29] GeV/c2. 

In less than 2% of the events, multiple B'^ candidates 
are present, and in these cases we choose that with a re- 
constructed D mass closest to the nominal mass value [9] . 
If more than one 5+ candidate share the same D can- 
didate, we select that with the smallest |Ai?|. In the 
following we refer to the selected candidate as Bsig- AH 
charged and neutral reconstructed particles not associ- 
ated with -Bsig, but with the other B decay in the event, 
-Bother; are called the rest of the event. 



III. BACKGROUND CHARACTERIZATION 

After applying the selection criteria described above, 
the remaining background is composed of non-signal BB 
events and continuum events. Continuum background 



events, in contrast to BB events, are characterized by a 
jet-like topology. This difference can be exploited to dis- 
criminate between the two categories of events by means 
of a Fisher discriminant F , which is a linear combination 
of six variables. The coefficients of the linear combina- 
tion are chosen to maximize the separation between sig- 
nal and continuum background so that T peaks at 1 for 
signal and at —1 for continuum background. They are 
determined with samples of simulated signal and contin- 
uum events, and validated using off-resonance data. In 
the Fisher discriminant we use the absolute value of the 
cosine of the angle between i?sig and i?othor thrust axes, 
where the thrust axis is defined as the direction maxi- 
mizing the sum of the longitudinal momenta of all the 
particles. Other variables included in F are the event 
shape moments Lq = SjPii ^^'^ ^2 = Z^iKI cos6'i|2, 
where the index i runs over all tracks and energy de- 
posits in the rest of the event; pi is the momentum; and 
9i is the angle with respect to the thrust axis of the i?sig- 
These three variables are calculated in the CM system. 
We also use the distance between the decay vertices of 
Ssig and D, the distance of closest approach between 
K meson tracks belonging to signal decay chain, and 
I At I, the absolute value of the proper time interval be- 
tween the i?sig and Bother decays [22]. The latter is cal- 
culated using the measured separation along the beam 
direction between the decay points of Bgig and -Bother 
and the Lorentz boost of the CM frame. The -Bother de- 
cay point is obtained from tracks that do not belong to 
the reconstructed -Bsigj with constraints from the -Bsig 
momentum and the beam-spot location. We use toes 
and F to define two regions: the fit region, defined as 
5.20 < toes < 5.29GeV/c2 and -5 < J" < 5, and the 
signal region, defined as 5.27 < toes < 5.29 GcV/c^ and 
< J" < 5. 

The BB background is divided into two components: 
non-peaking (combinatorial) and peaking. The latter 
consists of B-meson decays that have a well-pronounced 
peak in the toes signal region. One of the decay chan- 
nels which can mimic opposite-sign signal events, is the 
£?+ — > decay with D — > K+K^ and /9+ — > 7r+7r°. 

In order to reduce this contribution, we veto events 
for which the invariant K'^K~ pair mass ttix+k- is 
|m^+;^'- — A/o(PDQ) I > 20 MeV/c^ (with the D meson in- 
variant mass, Mb(vdg): taken to be 1864.83 MeV/c^ [9]). 
Simulations indicate that the remaining background is 
negligible. 

Another possible source of peaking BB background is 
the decay B^ — >■ Dn^ with D K~^Tr~Tr'^, which can 
contribute to the signal region of the same-sign sample 
due to the misidentification of the 7r+ as a AT"*". The 
number of events is expected to be about 8% of the total 
same-sign signal sample (see Table I). 

The charmless B~^ — >■ AT+AT^Tr+Tr" decay can also con- 
tribute to the signal region. The branching fraction of 
this decay has not been measured. Therefore the size of 
this background is estimated from the sidebands of the 
reconstructed D mass, 1.904 < Md < 2.000 GcV/c^ or 
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1.700 < Md < 1.824 GcV/c2. The result of the study is 
reported in Table I. In the final fit, we fix the yield 
of the same-sign BB peaking background to the sum 
of charmless and open-charm events. The opposite-sign 
background in the final event sample is assumed to be 
negligible. 

The overall reconstruction efficiency for signal events 
is (9.6 ± 0.1)% for opposite-sign signal events and (9.5 ± 
0.1)% for same-sign signal events. These numbers are 
equal within the uncertainty as expected. The composi- 
tion of the final sample is shown in Table I. 



IV. FIT PROCEDURE AND RESULTS 

To measure the ratios i?+ and R~ we perform extended 
maximum-likelihood fits to the toes a-nd J- distributions, 
separately for the B~^ and data samples. We write 
the extended likelihood functions as: 

with /±(x 1 9,N') = ^ ( ''"r;;^r°' 4,o.(x|g.ig,os)+ 

%ri/|g,ss(x|esig,ss)+E.<'/|.(x|e)); 

where /sig,ss(x| ^sig,ss), /sig,os(x |^'sig,os), and fB,{x.\d) are 
the probability density functions (PDFs) of the hypothe- 
ses that the event is a same-sign signal, opposite-sign 
signal, or a background event {Bi are the different back- 
ground categories used in the fit), respectively, N is 
the number of events in the selected sample, and N' 
is the expectation value for the total number of events. 
The symbol 6 indicates the set of parameters to be fit- 
ted. -^s±, total is the total number of signal events, 

= for the decays of the B^ meson, and 

is the total number of events of each background compo- 
nent. For the opposite-sign events the background comes 
from continuum and BB events. The peaking BB back- 
ground is introduced as a separate component in the fit 
to the same-sign sample. The fit is performed to the B^ 
sample (consisting of 15706 events) to determine i?+ and 
to the B~ sample (consisting of 15057 events) to deter- 
mine R~ . The PDFs for i?+ and R~ fits are identical. 
The i?ADS ratio is fitted to the same likelihood ansatz, 
but to the combined B'^ and B^ data sample. 

Since the correlations among the variables are neg- 
ligible, we write the PDFs as products of the one- 
dimensional distributions of toes and J-. The absence of 
correlation between these distributions is checked using 
MC samples. The signal toes distributions is modeled 
with the same asymmetric Gaussian function for both 
same-sign and opposite-sign events, while the distribu- 
tion is taken as a sum of two Gaussians. The continuum 
background toes distributions for the same and opposite- 
sign events are modeled with two different threshold AR- 



GUS functions [23] defined as follows: 



A{x)=x^l- (^-^^ •e'(' ), (10) 

where xq represents the maximum allowed value for the 
variable x, and c determines the shape of the distribu- 
tion. The Toes distribution of the non-peaking BB back- 
ground components are modeled with Crystal Ball (CB) 
functions that are different for same-sign and opposite- 
sign events [24] . The CB function is a Gaussian modified 
to include a power-law tail on the low side of the peak. 
The J- distributions for the BB background are approx- 
imated with sums of two asymmetric Gaussians. For the 
peaking BB background we conservatively use the same 
parameter set as for the signal. 

The PDF parameters are derived from data when pos- 
sible. The parameters for continuum events are deter- 
mined from the off-resonance data sample. The pa- 
rameters for the Toes distribution of signal events are 
extracted from the sample of i?+ — > Dtt^ with D 
ii'+7r~7r°, while for the parameters of the signal Fisher 
PDF we use the MC sample. The parameters of non- 
peaking BB distributions are determined from the MC 
sample. 

From each fit, we extract the ratios , or 

RadSi the total number of signal events in the sample 
(^B±,tot) along with the non-peaking background yields 
and threshold function slope for the continuum back- 
ground. Wc fix the number of peaking BB background 
events. 

To test the fitting procedure we generated 10000 
pseudo-experiments based on the PDFs described above. 
The fitting procedure is then tested on these samples. We 
find no bias in the number of fitted events for any com- 
ponent of the fit. Tests of the fit procedure performed on 
the full MC samples give values for the yields compatible 
with those expected. 

The main results of the fit to the data are summarized 
in Table II. 

The fits to the toes for > 0.5 and the J" distribution 
with Toes > 5.27 GeV/c^ are shown in Fig. 2, for the com- 
bined i?+ and B~ sample. These restrictions reduce the 
background and retain most of the signal events. Fig. 3 
shows the fits for the separate B~^ and B~ samples. 



V. SYSTEMATIC UNCERTAINTIES 

We consider various sources of systematic uncertain- 
ties, listed in Table III. One of the largest contributions 
comes from the uncertainties on the PDF parameters. 
To evaluate the contributions related to the toes and T 
PDFs, we repeat the fit varying the PDF parameters for 
each fit species within their statistical errors, taking into 
account correlations among the parameters (labeled as 
"PDF error" in Table III). 
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TABLE I: Composition of the final selected sample as evaluated from the MC samples normalized to data and from data for 
the charmless peaking background. The signal contribution is estimated using values of branching fractions from the PDG [9] 
and tb = 0.1 [10]. The errors are from the statistics of the control samples only. 



Sample 



Region Signal BB non-peaking Continuum Dtv Charmless peaking 



Same sign Fit 2252 ± 20 459 ± 12 

Signal 1921 ±18 147 ± 8 

Opposite sign Fit 28.7 ± 0.2 434 ± 12 

Signal 24.4 ±0.2 65 ± 5 



7403 ± 62 176 ± 14 28 ± 14 

203 ± 10 130 ±14 21 ± 14 

21201 ±104 - -2 ±9 

612 ±18 - -2 ±9 



I 20 
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- Full 
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BB peak 




BB non-peak 




~ continuum 





" Fisher 



Fisher 



FIG. 2: (color online) Distribution of (a,b) mES (with T > 0.5) and (c,d) T (with jties > 5.27 GeV/c^) and the results of 
the maximum likelihood fits for the combined _B"*" and B~ samples (extracting Rads), for (a,c) opposite-sign and (b,d) same- 
sign decays. The data are well described by the overall fit result (solid blue line) which is the sum of the signal, continuum, 
non-peaking, and peaking BB backgrounds. 
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FIG. 3: (color online) Projections of the 2D likelihood for mES with the additional requirement T > 0.5, obtained from the fit 
to the B^ (left) and B~ (right) data sample for opposite-sign events (extracting 7?^ and R~). The labeling of the curves is 
the same as in Fig. 2. 



To evaluate the uncertainties arising from peaking 
background contributions, we repeat the fit varying the 
the peaking BB background contribution vifithin its sta- 
tistical uncertainties and the errors of branching frac- 
tions, B, used to estimate the contribution. For the 
opposite-sign events only the positive part of the proba- 
bility distribution is used in the evaluation. 

Differences between data and MC (labeled as "Simu- 
lation" in Table III) in the shape of the T distribution 
are studied for signal components using the data control 
samples of _B+ — )• Dn^ with D A'+7r~7r". These pa- 



rameters are expected to be slightly different between the 
B — Dtt and B — DK samples. We conservatively take 
the systematic uncertainty as the difference in the fit re- 
sults from the nominal parameters set (using MC events) 
and the parameters set obtained using the B — > Dtt data 
sample. 

The systematic uncertainty attributed to the cross- 
feed between opposite-sign and same-sign events has been 
evaluated from the MC samples. The number of same- 
sign events passing the selection of the opposite-sign 
events is taken as a systematic uncertainty. The efficien- 
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TABLE II: Results of fits to the B+, B", and the combined 
and B~ samples, including the extracted number of signal 
and background events and their statistical errors. 



Sample 


B+ 


B- 


B+ and B~ 


R, 10"^ 


c + 12 
•-'-10 


,2 + 12 






1032 ± 41 


946 ± 39 


1981 ± 57 


/ybkg 
SB, OS 


305 ± 52 


120 ± 36 


402 ± 65 


/ybkg 


315 ± 44 


329 ± 44 


644 ± 62 


■'^cont.OS 


10290 ±111 


10017 ± 105 


20329 ± 154 


A^bkg 

-'^cont.SS 


3660 ± 69 


3539 ± 68 


7203 ± 76 



TABLE III: Systematic errors for R and Raus in units of 
10"^ 



Source 


R+ 


R- 


-Rads 


PDF error 


+1.1 

-1.8 


1.1 


1.0 


Same sign peaking background 


0.2 


0.5 


0.2 


Opposite sign peaking background 


+ 
-3.C 


+0 
-3.6 


+0 
-3.4 


Simulation 


0.6 


0.6 


0.7 


B errors 


0.2 


0.6 


0.4 


Crossfeed contribution 


0.1 


0.4 


0.3 


Efficiency ratio 


0.1 


0.4 


0.3 


Combined uncertainty 


+ 1.2 

-4.1 


+ 1.6 
-3.9 


+1.4 
-3.7 



cies for same-sign and opposite-sign events were verified 
to be the same within a precision of 3% [25]. We hence 
assign a systematic uncertainty on based on varia- 
tions due to changes in the efficiency ratio by ±3%. 

The systematic uncertainties for the ratios R^, R~, 
and i?ADS are summarized in Table III. The overall sys- 
tematic errors represent the sum in quadrature of the 
individual uncertainties. 



VI. EXTRACTION OF re 

Following a Bayesian approach [26], the probability 
distributions for the i?+ and ratios obtained in the 
fit are translated into a probability distribution for rs 
using Eqs. 8 and 9 simultaneously. We assume the fol- 
lowing prior probability distributions: for ru a Gaussian 
with mean 4.7 x I0~^ and standard deviation 3 x 10""^ [9] ; 
for and we use the likelihood obtained in Ref. [8], 
taking into account a 180 degree difference in the phase 
convention for do ; for 7 and S b we assume a uniform dis- 
tribution between and 360 degrees, while for a uni- 
form distribution in the range [0, 1] is used. We obtain 
the posterior probability distribution shown in Fig. 4. 
Since the measurements are not statistically significant, 
we integrate over the positive portion of that distribution 
and obtain the upper limit tb < 0.13 at 90% probability. 



and the range 

TB e [0.01,0.11] at 68% probability, (11) 
and 0.078 as the most probable value. 



xlO'^ 




0.1 0.2 0.3 0.4 0.5 



FIG. 4: Bayesian posterior probability density function for 
rs from our measurement of R'^ and R^ and the hadronic 
D decay parameters ro, Sd, and ko taken from [8] and [9]. 
The dark and light shaded zones represent the 68% and 90% 
probability regions, respectively. 



VII. SUMMARY 

We have presented a study of the decays D^K^ 
and D°K^, in which the D° and Z)° mesons decay 

to the A'^tt^tt" final state using the ADS method. The 
analysis is performed using 474 x 10^ BB pairs, the full 
BABMi dataset. Previous results [7] are improved and 
superseded by improved event reconstruction algorithms 
and analysis strategies employed on a larger data sample. 

The final results are: 

= (5tl^(stat):^:^(syst)) x 10-^, (12) 
R- = (^12ijg(stat)+4(syst)) x 10~^ (13) 
i?ADS = (9.lt?j(stat)+J:^(syst)) x (14) 

from which we obtain 90% probability limits: 

R+ < 23 X 10-3, (15) 
R~ < 29 X 10-^ (16) 
i?ADS < 21 X 10-3. (17) 

From our measurements we derive the limit 

rB < 0.13 at 90% probability. (18) 
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